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Abstract. It is widely believed that gamma-ray bursts are produced by a jet-like outflows directed towards the 
observer, and the jet opening angle(0j) is often inferred from the time at which there is a break in the afterglow 
light curves. Here we calculate the GRB afterglow light curves from a relativistic jet as seen by observers at a 
wide range of viewing angles (6'„) from the jet axis, and the jet is uniform or non-uniform(the energy per unit 
solid angle decreases smoothly away from the axis e{e) <x {6/6c)~''). We find that, for uniform jet(fc = 0), the 
afterglow light curves for different viewing angles are somewhat different: in general, there are two breaks in the 
light curve, the first one corresponds to the time at which 7 ~ (Sj — and the second one corresponds to 

the time when J ^ [Oj + 9^)~^ . However, for non-uniform jet, the things become more complicated. For the case 
d-u = 0, we can obtain the analytical results, for k < 8/{p + 4) (where p is the spectral index of electron energy 
distribution) there should be two breaks in the light curve correspond to 7 ~ 6^"^ and 7 ~ respectively, while 
for k > 8/{p + 4) there should be only one break corresponds to 7 ~ and this provides a possible explanation 
for some rapidly fading afterglows whose light curves have no breaks since the time at which 7 ~ is much 
earlier than our first observation time. For the case 6^ 7^ 0, our numerical results show that, the afterglow light 
curves are strongly affected by the values of 9v, 9c and k. If 0„ is close to 9c and k is small, then the light curve is 
similar to the case of fe = 0, except the fiux is somewhat lower. However, if the values of dv/9c and k are larger, 
there will be a prominent flattening in the afterglow light curve, which is quite different from the uniform jet, and 
after the flattening a very sharp break will be occurred at the time 7 ~ (6v + 9c)~^. 
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1. Introduction 

Gamma-ray bursts (GRBs) are known as an explosive 
phenomenon occurring at cosmological distances, emitting 
> large amount of energy mostly in the gamma-ray range 
■ (see, e.g. Piran 1999; Cheng & Lu 2001 for a review). 
Observations show that some of GRBs are emitting an 
extremely large energy with 3> 10^^ ergs if emission 
is isotropic. For example, GRB990123, the most energetic 
GRB event detected so far, has an isotropic gamma-ray 
energy of ^ 3.4 x 10^'* ergs, which corresponds to the rest- 
mass energy of ^ I.9M0 (Kulkarni et al. 1999). Such a 
crisis of extreme large energy forced some people to think 
that the GRB emission must be highly coUimated in order 
to reduce the total energy. 

The second reason for GRB emission being jet-like 
comes from the fact that there are sharp breaks in the light 
curves of some GRBs' afterglows, such as GRB990123 
(Kulkarni et al. 1999; Castro-Tirado et al. 1999) and 
GRB990510 (Harrison et al. 1999; Stanek et al. 1999), etc.. 
These observed breaks have generally been interpreted as 
evidence for coUimation of the GRB ejecta, since Rhoads 
(1999) and Sari et al. (1999) have pointed out that the lat- 



eral expansion of the relativistic jet can produce a sharp 
break in the afterglow light curve. 

However, in the current afterglow jet models, it is gen- 
erally assumed that the jet is uniform, and the line-of-sight 
is just along the jet axis. It is obvious that these assump- 
tions are usually not true, since some GRB models predict 
that the jet may be non-uniform, within which the energy 
per unit solid angle decreases away from the jet axis, such 
as e oc 6^'^ (e.g. MacFadyen et al. 2001), and the proba- 
bility that the line-of-sight just crosses the jet axis is near 
zero. So it is very important to investigate the situation 
that the jet is non-uniform and is seen by observers at a 
wide range of viewing angles from the jet axis. Some au- 
thors have already considered the situation of anisotropic 
jets (Meszaros, Rees & Wijers 1998; Salmonson 2001; Dai 
& Gou 2001; Zhang & Meszaros 2002; Rossi, Lazzati & 
Rees 2002). In this paper, we give a detailed calculation 
of the emission from anisotropic jets, including the effect 
of equal-arrival-time surface. In the next section we dis- 
cuss the dynamical evolution of the jet, in section 3 we 
calculate the emission from uniform jet for different view- 
ing angle, in section 4 we calculate the emission features 
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from non-uniform jet, and finally we present some discus- 
sions and conclusions. 



2. Dynamical evolution of the jet 

Now wc consider an adiabatic relativistic jet expanding 
in the surrounding medium. For energy conservation, the 
evolution equation is 



(1) 



where 7 is the bulk Lorentz factor, A A'' = ^nr^AQ is the 
particle numbers swept by the jet within a solid angle Afi, 
n is the surrounding medium density, and e is the energy 
per unit solid angle. 

It is well known that for relativistic blast waves, the 
received photons at time T are not emitted at the same 
time. A photon that is located at radius r and is emitted 
with an angle from the linc-of-sight will reach the ob- 
server at time T = r(l — /3/Lt) / (/3c), where fi = cos Q. Then 
we can obtain the jet evolution 



7 = (^)l/2(^pc2)-l/2(_Ml^-3/2 



(2) 



We sec that the jet Lorentz factor evolution is dependent 
on the angle Q, for different values of 6 the relation be- 
tween 7 and T are different. Here we neglect the sideways 
expansion of the jet since this process is very unclear. For 
relativistic jet, 



7(1-/?/.)-^/^- 3.4 xlO^(^) V^T-„f 



(3) 



where 654 is the energy in units of 10^^ ergs, ni is the 
surrounding density in units of 1 atom cm^^, and T^ay is 
the observed time in units of 1 day. In the case of 7 ^ 1 
and 6 <i; 1, we have 



7(i + 7'e^) 



2cv2^-3/8 



V > day 



(4) 



This equation can be solved numerically, and by fitting 
the numerical results, we obtain the analytical solution 



_ / ioi-2^(ff )i/«r-„f , for 7e < 1 



7 = 



10 



5.08 



m f'^T,^T^\ for 76 > 1 



(5) 



It is obvious that the usual solution 7 oc r-3/8 is valid 
only when 8 < 7"^, while when Q > 7"^ the relation 
becomes 7 oc T~^/^. 

3. emission from uniform jet 

Emission features from uniform jet has been discussed by 
many authors, and it is widely believed that there is a 
sharp break in the GRB afterglow light curve corresponds 
to the time 7 ~ 9j^, where 9j is the jet half-opening angle. 
However, this is true only when the observer line-of-sight 
just crosses the jet axis, and in fact this probability is very 
small. So here we calculate the jet emission for various 
viewing angles. 



We follow our previous paper (Wei & Lu 2000) to cal- 
culate the emission flux from the jet. We assume that 
the line-of-sight is along z-axis, the symmetry axis of 
the jet is in the y — z plane, 9y is the angle between 
the line-of-sight and the symmetry axis, and the radi- 
ation is isotropic in the comoving frame of the ejecta. 
In order to see more clearly, let us establish an auxil- 
iary coordinate system {x',y',z') with the 2;'-axis along 
the symmetry axis of the cone and the x' parallel the x- 
axis. Then the position within the cone is specified by 
its angular spherical coordinates and (f) < 9 < 6j, 
< 4> < 2tt). It can be shown that the angle O between a 
direction {9, <p) within the cone, and the line-of-sight satis- 
fies cos6 = cos9yCos9 — sin9ySin9sin(f). Then the observed 
flux is 



sin9d9D^l'{uD-'^] 



cP 



(6) 



where D = [7(1 — /3cos6)]~^ is the Doppler factor, (} = 
(1— 7~^)^/^, V = Du', I'{v') is the speciflc intensity of syn- 
chrotron radiation at v' , and d is the distance of the burst 
source. Here the quantities with primes are measured in 
the comoving frame. 

It is generally believed that the electrons have been 
accelerated by the shock to a power law distribution 
ne(7e) oc 7"^" for 7„,„ < 7e < ^max, and consider the 
synchrotron radiation of these electrons, we can obtain 
the observed flux 



o p — 1 



2tt pSj 

d(l) / sin9d9f{j){l- (3^)- 
Jo 



(7) 



where T is the time measured in the observer frame 

fh) = i)(^+^^/'2/(7)(^+')/SLn/3'> y(7) = 

ji, 7 is just the ratio of specific heats, "imAn = 6; (7 ~ 
P p-2 

e p-1 

relativistic case, we have 



77+1 
7 



1)^|_^ mp(me) is the mass of proton (electron). In the 



oc T^z/- — 



d(j) I sin9d9y 



.p-3 



{l-Pn)-"^ (8) 



Combining the jet evolution results (eq. (2)), finally we 
obtain the observed flux 



F.ocT-tf+^t.-^ 



27r 



sin9d9e^{l-l3fj,)P-'^°{9) 



This is the basic equation for our calculation. Here we take 
9j = 0.1, p — 2.5 and ni = 1 throughout this paper. 

For uniform jet, the energy per unit solid angle e is 
independent of 9, cf), we take €54 = 1. Then using equations 
(2) and (9), we can calculate the jet emission for different 
viewing angles. Fig.l gives our numerical results. From 
Fig.l we see that the afterglow light curves for different 
viewing angles are somewhat different: in general, there 
are two breaks in the light curve, the flrst one corresponds 
to the time at which 7 ^ {9j — 9v)^^, and the second one 
corresponds to the time when 7 ~ {9j + 9y)^^ . This is 
quite different from the previous results, which think there 
is one break occurred at 7 ~ 
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Fig. 1. The afterglow light curves from uniform jet for various Fig. 2. The afterglow light curves from non-uniform jet for the 
viewing angles. case 9^ = 0. 



4. Emission from non-uniform jet 

Uniform jet is only a special case, in general we expect 
the jet is non-uniform, for instance, the coUapsar model 
predicted that the jet is non-uniform, within which the 
energy per unit solid angle decreases away from the jet 
axis (e.g. MacFadyen et al. 2001). Here we suppose the 
energy per unit solid angle has the form 



£0 



for 
for 



< 6* < 6 



(10) 



where 6c is introduced to avoid a divergence at — Q. For 
arbitrary viewing angle 9y, the situation is very compli- 
cated, so first we consider the case ^t, = 0, since this can 
give an analytical analysis. 

In the case of 6'i, = 0, we have /i = cos 9. It should 
be noted that there is an interesting phenomenon, the 
values of ^9 increase with the angle 0, i.e. there is a 
characteristic angle 9i,, at which the value 7*0* = 1, 
when 9 < 9i, the value 76* < 1, and when 9 > 9i, 
the value ^9 > 1. It is obvious that the main con- 
tribution of emission comes from the region 9 < 9^,, 
so 9i, is an important quantity. It can be shown that 
when T < Ti = l-7(^)^/^(o%)^/^ hours, the value 
9. = 10-i-27(£Mi)-i/8T3/8 < e^, and when T < T2 = 
5.2 X (0.2)'=/3(lMl)i/3(_|_)fe/3( e^)(8-fe)/3 days, the valuc 

Q _ 1 Q-1.27x8/(8-fc)[' eo.54 \-l/(8-fc)fl-fc/(8-fc) ji3/(8-fc) ^ 
* V ni '' ^ day — 

Taking equation (5) and the approximate expression 
1 - ^/.t ~ (272)-! for 9 < 9^, or 1 - ~ ^9^ for 9 > 9^, 
we can get the analytical results: (1) when T < Ti, the 
observed flux oc T-^ip-^)/'^- (2) when Ti < T < T2, 

3(2p+fc-2) , 

the flux (X T ^ — for fc < or F^, cx T~3p/4 

for A: > (3) when T > T2, the flux F^ oc T-^p/^. 

From this result we see that, for smaller value of k, the 



to 



3£ 

4 



3(P-1) 
4 

is gradual and smooth with a timescale of about 



transition of the afterglow light curve index from ^ 



T2, while for larger values of k, the transition is rapid with 
a timescale of about Ti . 

In order to verify the above results, we also make a 
numerical calculation for the case 9v = 0. Fig. 2 gives our 
results. It is obvious that the numerical results are consis- 
tent with the analytical results, for larger value of fc, the 
steepening of the light curve is more rapidly. We suggest 
this may explain some afterglow light curves which decay 
rapidly and have no breaks, since for larger value of k the 
transition time (~ Ti) is earlier than our first observation 
time. 

However, it should be noted that the appearance of the 
early break in the light curve (corresponding to the time 
when J ^ 9^^) is due to the assumed energy distribution 
function (equation (10)), and the sharpness of this break is 
primarily dependent on the discontinuity in slope of e for 
the idealized model of equation (10) at = 9c- It is obvious 
that for a realistic energy distribution, the transition of e 
from roughly constant for 9 < 9c to e (x 9^^ for 9 > 9c 
should be smooth. Therefore this break may be washed 
out by a realistic energy distribution. 

For the general case 9^ ^ 0, we calculate the afterglow 
light curves numerically using equations (2), (9) and (10), 
the results are given in Fig. 3 - Fig. 6. From these figures we 
see that the afterglow light curves are dependent on the 
values of 9^, 9c and k. When 9^ is near 9c and k is small, 
then the light curve is similar to the case of 0„ =0, except 
the fiux is somewhat lower. However, if the viewing angle 
9y is larger than 9c and k is not very small, then there 
will be a prominent flattening in the afterglow light curve, 
which is quite different from the case of uniform jet, and 
after the flattening a very sharp break will be occurred 
around the time 7 {9^ +9c)~^ ■ If the viewing angle 9y is 
larger than the jet half-opening angle 9j, then the flux will 
first increase with time until the Lorentz factor is about 
1 ^ {Sv ~ 9j)^^ , thereafter the flux begin to decay with 
time. 
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Fig. 3. The afterglow light curves from non-uniform jet for the 



= 0.06. 
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Fig. 4. The afterglow light curves from non-uniform jet for the 
case 6v = 0.1. 

5. Discussion and conclusion 

In this paper we calculate the GRB afterglow light curves 
from relativistic jets in more details, assuming that the jet 
may be uniform or non-uniform, and the observer locate 
at arbitrary angle with respect to the jet axis. We have 
shown that there are several distinct features of our jet 
emission compared with previous jet model. 

In previous analysis, it is generally assumed that the 
jet is uniform and the line-of-sight is just along the jet 
axis, in this case the afterglow light curves have a break 
at the time 7^6*^^. However, if the viewing angle 9^ ^ 0, 
we have shown that there should be two breaks in the 
light curve, the first one corresponds to the time at which 
^ ^ {9j — 6y)^^ ^ and the second one corresponds to the 
time when 7 ~ {9j + 6'„)^^, although these transitions are 
very smoothly. 

If the jet is not uniform, within which the energy dis- 
tribution is given by eq.(lO), then the calculation is more 
complicated. In the case of ^i, = 0, we can give an analyti- 
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Fig. 5. The afterglow light curves from non-uniform jet for the 
case 6^ = 0.2. 
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Fig. 6. The afterglow light curves from non- uniform jet (k — 2) 
for different viewing angles. 

cal result, for k < 8/ (p+4) (where p is the spectral index of 
electron energy distribution) there should be two breaks 
in the light curve correspond to 7 ^ 9^^ and ^ ^ 9J^ 
respectively, while for fc > 8/ (p -t- 4) there should be only 
one break corresponds to 7 ^ 9^^, after that the flux 
decays as Fi, oc T^Sp/**^ \Ye argue that this may explain 
some rapidly fading afterglows whose light curves have no 
breaks, since the time Ti, at which 7 ~ 9~^, is usually 
earlier than our first observation time. 

If the jet is non-uniform and the viewing angle 9^ ^ 0, 
only numerical results can be given. We have shown that 
in this case the shape of afterglow light curve is dependent 
on the values of 9^, 9c and k. When 9^ is near 9c and k is 
small, then the light curve is similar to the case of 9^ = 0, 
except the fiux is somewhat lower. However, if the viewing 
angle 9y is larger than 9c and k is not small, then there 
will be a prominent flattening in the afterglow light curve, 
which is quite different from the case of uniform jet, and 
after the flattening a very sharp break will be occurred 
around the time j ^ {9v + 9c)~^ ■ We think this is a main 
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difference between the uniform and non-uniform jet, and 
we can identify whether the jet is uniform or not by this 
feature. 

It is not very difficult to understand why sometimes 
there is a flattening in the Ught curve. It is weU known 
that, for a rclativistic blast wave with Lorcntz factor T ^ 
1, the observer can only observe a solid angle around 9y 
with a half opening angle of order T~^, the contribution 
from other components can be neglected, so when the blast 
wave decelerates, the observer can see larger components. 
For the case of non-uniform jet and 6y > 6c, the energy 
at 9 = 9y is much smaller than that of smaller 6, so when 
r decreases, the observer can observe more region with 
larger energy (smaller 6), so there will be a flattening in 
the light curve. 

We would point out that, in our calculation we have 
neglected the sideways expansion of the jet since this pro- 
cess is too complicated. However, wc know that, in fact 
this process is an important issue in determining the shape 
of a light curve, since it will likely significantly change the 
shape of the light curve for both the uniform and non- 
uniform jets. So we suggest that for a more realistic cal- 
culation the sideways expansion should be taken into ac- 
count. 
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